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Pion form factor in QCD sum rules, local duality approach, and 0(^2) 
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Using the results on the electromagnetic pion Form Factor (FF) obtained in the 0(aa) QCD sum rules with 
non-local condensates [T] we determine the effective continuum threshold for the local duality approach. Then we 
apply it to construct the 0(0^) estimation of the pion FF in the framework of the fractional analytic perturbation 
theory. 
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In this paper we use the results on the electro- 
magnetic pion FF F-^iCf'') in the spacelike region 
Q2 = 1 _ 10 GeV^, obtained in the 0{as) QCD 
Sum Rule (SR) approach with nonlocal conden- 
sates (NLC) [I], in order to estimate the next-to- 
next-to-leading (or the 0{o?^)) order (NNLO) cor- 
rection to the pion FF using both the Fractional 
Analytic Perturbation Theory (FAPT) and the 
Local Duality (LD) approach. First, we describe 
the LD approximation in the pion FF calculation 
and discuss how it is possible to obtain the main 
ingredient of this approach, namely, the contin- 
uum threshold so(Q^) as a function of . Then 
we explain why one need to use FAPT in esti- 
mating the NNLO result for the factorized part 
of the pion FF. Next, we describe our model [5] 
for matching function, gluing soft and hard parts 
to the complete pion FF. After that we show how 
to apply FAPT to estimate the NNLO pion FF 
using improved model for matching function. 



1. Pion FF in the Local Duality approach 



The LD SR [3|4| is produced from the original 
QCD SR in the oo limit. For this reason it 

has no condensate contributions. The main non- 
perturbative ingredient in this approach is the ef- 
fective continuum threshold Sq^ — it inherits all 
the nonperturbative information from the origi- 



nal QCD SR. At the {I + l)-loop order we have 
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where S should be substituted by the LD effective 
threshold, s]^'^^\Q'^), and p^^\si, S2,Q'^) is the 
three-point [l -I- l)-loop spectral density. In the 
leading order the integration can be done analyt- 
ically: F^^-'^''\q \ S) = [S/{A^ ^m] [1 - (Q^ + 
65) /(Q2 + 45) V<3V(Q' + 45)]. The LD pre- 
scription for the corresponding correlator ^5 4^ im- 
plies the relations 



47r2/x - 0-7 GeV^ 



and 



LD;(i; 



(0) 



l-|-a,(g2)/7r 



0.6 GeV^ 



where Qq is of the order of Sq^'^°''(0). This pre- 
scription is a strict consequence of the Ward iden- 
tity for the AAV correlator due to the vector- 
current conservation. In principle, the depen- 
dence of the LD parameter s}^^{Q^) ([1]) should 
be determined from the QCD SR at > 
1 GeV^. But as explained in [1] the standard 
QCD SR becomes unstable at > 3 QgV^ 
because of the appearance of terms in the con- 
densate contributions linearly growing with 
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[6l7j . For this reason, this dependence was 
known only for < 3 GeV^ and, therefore, 
most authors usually used the constant approx- 



imation s 



LD;(0) 
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^0 



(0), like in [H 



or a slightly Q -dependent approximation 

s^^'^W') =^ 4 7r2/^/(l + as(g2)/7r), likein [10]. 

But now, due to the knowledge of the NLC 
QCD SR prediction [1^ for the pion FF for ^ 
1 — 10 GeV^, we can estimate the effective LD 
threshold Sq^{Q^) which reproduce these predic- 
tions in the LD approach. Results are shown in 
Fig. 1. It can be represented in this Q"^ range by 
the following interpolation formula: 



s^^(Q-^) = 0.57 + 0.461 tanh 



0.0954 Q 

GeV^ 



21 
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We see that Sq^ {Q^ ) in the mentioned range of 
is monotonically increasing function. Therefore 
sLD(q2) ^ sSR(Q2) ^ 0.7 GeV2 and due to this 
difference the LD approaches of [ 219110] produces 
significantly lower predictions for (5^i^7r(Q^) as 
compared with QCD SRs with NLC. 
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Ft gure 1. The effective LD threshold as a 
function of Q^. 



2. Using FAPT for NNLO estimation 

In order to estimate the NNLO contribution 
to the pion FF in the QCD SR approach one 
needs to know the three-loop spectral density 
pI^^ (si, S2, Q^) and this is a complicated task. We 
want to avoid this tricky calculation and suggest 
to use the known coUinear two- loop result, the LD 



model for the soft part with improved Sq^{Q^), 
and matching procedure of 2 . And one needs to 
use (F)APT for the two-loop coUinear expression, 
as been shown in [2)llj . in order to have practi- 
cal independence with respect to renormalization 
and factorization scale setting. 

To 'glue' together the LD model for the soft 
part, F^^''"^\q'^) (which is dominant at small 
< 1 GeV^), with the perturbative hard- 
rescattering part, F^'^'^^''^^\q^) (which provides 
the leading perturbative 0{as) + 0{al) correc- 
tions and is dominant at large 3> 1 GeV^), 
in such a way as to ensure the validity of the 
Ward identity (WI) F^^'^^\o) = 1, we apply the 
matching procedure, introduced in [5]: 



FWI;(2)(Q2)^j^LD,(0)(q2 



2s, 



(2) 



^pQCD,(2)(q2) (4) 



with Sq ~ 0.6 GeV^. In order to test the qual- 
ity of the matching prescription given by Eq. ([J]), 
we propose to compare it with the LD model ([T]) 
evaluated at the one-loop order (i.e., in the 0{as)- 
approximation |9ll0j ). To this end, we construct 
the analogous 0(as)-model F^^''^^\q'^), where 
we substitute FpQ^°'(')(Q2) by F^'^^'^'^^^^Q^) = 

2as((5^) SQ^'^°\o)/Tr Q'^ and imply the same pre- 
scription for the effective LD threshold as in [TU] , 
i. e. (l2bD. It is worth to note that the model 



^wi,(i) (^q2-j^ suggested in 2 , works quite well, al- 
though it was proposed without the knowledge of 
the exact two-loop spectral density, which became 
available later [S] . The key feature of this match- 
ing recipe is that it uses information on Ft^[Q^) 
in two asymptotic regions: 

1. 0, where the Ward identity dic- 
tates F^(0) = 1 and hence Fj^iQ^) ~ 

f^'^'^'Hq^ 

2. Q2 ^ ^iiere F^{Q^) ~ F^^™'^^\q^) 

in order to join properly the hard tail of the 
pion FF with its soft part. Numerical analy- 
sis shows that the applied prescription yields a 
pretty accurate result, with a relative error vary- 
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ing in the range 5% at ^ ^ ggv^ to 9% at 
q2 = 3_3o GeV2. 

Now, when the spectral density P3^^(si, S2, Q^) 
is known [9j, it is possible to improve the repre- 
sentation of the LD part by taking into account 
the leading 0{as) correction in the electromag- 
netic vertex. We suggest the following improved 
WI model 

fZ'£\q^s) = f^^-^(^HQ',s) 

S as{Q^) ( 2S y 

with subsequent substitution S — > s^^'^^\q'^). 
We explicitly display the dependence on the 
threshold S in Eq. ©—the aim being to apply 
it later on with S = Sq^(Q^), the latter value 
being extracted by comparing the NLC QCD SR 
results with the LD approximation. Numerical 
evaluation of this new WI model in comparison 
with the exact LD result in the one-loop approx- 
imation shows that the quality of the matching 
condition is improved: the relative error is re- 
duced, reaching 4% at = 1 _ 10 GeV^. 

Proceeding along similar lines of reason- 
ing, we construct the two-loop WI model 

tain 

, S as{Q') ( 2S y 

+ ^ FPAPT,(2)/02n ( Ql \ /gN 

where F^^^'^'''^\q'^) is the analyticized expres- 
sion generated from i^^*^^^'*-^-* (Q^) using FAPT 
(see Refs. [12113111) ) to get a result which appears 
to be very close to the outcome of the default scale 
setting (/ip^ = ^1 = Q2-)^ investigated in detail in 
[2] in the APT approach. FAPT is needed here in 
order to obtain analytic expressions for the pion 
FF in both possible cases of factorization scale 
setting: 



(i) For fip — Q2 there appear factors of the 
type [q!s(Q^)] with fractional powers ly = 
7„/(2 6o) due to the pion distribution am- 
plitude evolution; 

(ii) For /i| = const there appears factor 

In any case, the NNLO correction involves the 
analytic image of the second power of the cou- 
pling, -42 (Q^). For this reason we name the whole 
^FAPT,(2) ^q2-j term as 0{A2) contribution. 
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Figure 2. We show as a narrow dash-dotted 
strip the predictions for the pion FF, obtained 
in the two- loop WI model, Eq. ([6]), using the im- 
proved Gaussian model. The width of the strip 
is due to the variation of the Gegenbauer coeffi- 
cients 02 and 04 (needed to calculate the coUinear 
part F^'^'^^''"^\q'^)) in the corresponding shaded 
bands for the pion DA (indicated by the central 
solid line). Note that this dash-dotted strip shows 
the effect of 0(^2) correction only for the central 
solid curve of the shaded band. 

Interesting to note here, that in the case of the 
one-loop approximation the relative error of WI 
model (O appears to be of the order of 10%. The 
relative 0(as)-contribution to the pion FF is of 
the order of 10%, as has been shown in [2111] . 
Hence, the relative error of our estimate is of the 
order of 1% — provided we take into account the 
O(as)-correction exactly via the specific choice of 
soiQ^), as done in 
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The results obtained for the pion FF with our 
two- loop model, i.e., Eq. ([5]), and using the ef- 
fective LD thresholds Sq^{Q^), are displayed in 
Fig. 2. We see from this figure that the main 
effect of the NNLO correction peaks at > 
4 GeV^, reaching the level of 3 - 10%. 

Conclusions 

• We showed here that the local duality 
model for pion FF suffers from the thresh- 
old so(Q^) uncertainty. We fixed this un- 
certainty by demanding that the LD model 
with right setting for so(Q^) should repro- 
duce results for pion FF obtained in the 
Borel SRs with NLC [I]. Our results show 
that s^^(Q2) grows with Q^. 

• Our rough model for matching function 
in [2] appears to be of good quahty (w 10%) 
and we managed to improve it here to have 
quality of « 5%. 

• Using FAPT and improved matching func- 
tion we estimated NNLO correction to the 
pion FF to be of the order of w 3 - 10%. 

• Our strip of predictions for the pion FF is 
in a good agreement with existing experi- 
mental data of Cornell [T4l and JLab [15] . 

Acknowledgments 

This work was supported in part by the Russian 
Foundation for Fundamental Research, grants 
No. 07-02-91557, 08-01-00686, and 09-02-01149, 
the BRFBR-JINR Cooperation Programme, con- 
tract No. F08D-001, the Deutsche Forschungsge- 
meinschaft (Project DFG 436 RUS 113/881/0- 
1), and the Heisenberg-Landau Programme un- 
der grant 2009. 

REFERENCES 

1. A. P. Bakulev, A. V. Pimikov, and N. G. Ste- 
fanis, Phys. Rev. D79 (2009) 093010. 

2. A. P. Bakulev, K. Passek-Kumericki, W. 
Schroers, and N. G. Stefanis, Phys. Rev. D70 
(2004) 033014, 079906(E). 



3. V. A. Nesterenko and A. V. Radyushkin, 
Phys. Lett. B115 (1982) 410. 

4. A. V. Radyushkin, Acta Phys. Polon. B26 
(1995) 2067. 

5. M. A. Shifman, A. I. Vainshtein, and V. I. 
Zakharov, Nucl. Phys. B147 (1979) 385; ibid. 
448; ibid. 519. 

6. B. L. loffe and A. V. Smilga, Nucl. Phys. B216 
(1983) 373. 

7. V. A. Nesterenko and A. V. Radyushkin, 
JETP Lett. 39 (1984) 707. 

8. A. P. Bakulev, A. V. Radyushkin, and N. G. 
Stefanis, Phys. Rev. D62 (2000) 113001. 

9. V. V. Braguta and A. L Onishchenko, Phys. 
Lett. B591 (2004) 267. 

10. Victor Braguta, Wolfgang Lucha, and Dmitri 
Melikhov, Phys. Lett. B661 (2008) 354. 

11. A. P. Bakulev, Phys. Part. Nucl. 40 (2009) 
715. 

12. A. P. Bakulev, S. V. Mikhailov, and N. G. 
Stefanis, Phys. Rev. D72 (2005) 074014, 
119908(E); ibid. D75 (2007) 056005; ibid. D77 
(2008) 079901(E). 

13. A. P. Bakulev, A. I. Karanikas, and N. G. 
Stefanis, Phys. Rev. D72 (2005) 074015. 

14. C. J. Bebek et al, Phys. Rev. D9 (1974) 1229; 
ibid. D13 (1976) 25; ibid. D17 (1978) 1693. 

15. G. M. Huber et al, Phys. Rev. C78 (2008) 
045203. 



